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Abstract 

The s— wave meson-baryon scattering is analyzed for the strangeness S = — 1 and isospin 
1 = sector in a Bethe-Salpeter coupled channel formalism incorporating Chiral Symmetry. 
Four channels have been considered: ttT,, KN, r/A and AH. The required input to solve 
the Bethe-Salpeter equation is taken from lowest order Chiral Perturbation Theory in a 
relativistic formalism. There appear undetermined low energy constants, as a consequence 
of the renormalization of the amplitudes, which are obtained from fits to the ttT, — * 7rS 
mass-spectrum, to the elastic KN — ► KN and KN — > 7rS i-matrices and to the K~p — ► r/A 
cross section data. The position and residues of the complex poles in the second Riemann 
Sheet of the scattering amplitude determine masses, widths and branching ratios of the 
Sq\— A(1405) and — A(1670) resonances, in reasonable agreement with experiment. A good 
overall description of data, from 7r£ threshold up to 1.75 GeV, is achieved despite the fact 
that three-body channels have not been explicitly included. 
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1 Introduction 



Baryon resonances are outstanding features in elastic and inelastic meson-baryon scattering and 
signal the onset of non-perturbative physics. Constituent quark model approaches describe them 
as excited baryonic bound states, and the coupling to the continuum is obtained by evaluating 
transition matrix elements [1] but comparison with data can only be done once the scattering 
problem is solved. In such a scheme, the underlying quark constituent nature of hadrons is taken 
into account but implementation of Chiral Symmetry (CS) becomes difficult. In the region of 
low energies, it seems appropriate to start considering the hadrons as the relevant degrees of 
freedom where CS not only proves helpful to restrict the type of interactions between mesons 
and baryons, but also provides an indirect link to the underlying Quantum Chromo Dynamics 
(QCD) [2]. For processes involving Baryons and Mesons, Heavy Baryon Chiral Perturbation 
Theory (HBChPT) [3, 4] incorporates CS at low energies in a systematic way, and has provided 
a satisfactory description of nN scattering in the region around threshold [5, 6, 7]. It suffers, 
however, from known limitations. Firstly, the expansion is manifestly not relativistically invariant, 
and some convergence problems, specially for the scalar form factor, have been pointed out and 
solved by defining a suitable regularization scheme, the so-called infrared regularization (IR) [8]. 
The IR scheme has also successfully been applied to tcN elastic scattering [9]. The previous 
remarks concern a theory with only n and iV present. The A degrees of freedom have been 
explicitly included [10] within the HBChPT scheme, where the nucleon-A splitting is considered 
to be of order of the pion mass, pointing toward a better convergence for the nN scattering data 
than the case without explicit A's. More recently, the work of Ref. [11] has compared the HBChPT 
scheme to the IR one in the presence of explicit A's degrees of freedom for nN scattering, showing 
that HBChPT describes data up to much higher pion CM kinetic energies than the IR method. 

A second limitation of any of the previous approaches, both in the HBChPT as well as in 
the IR frameworks and with or without explicit A, is that they are based on a perturbative 
expansion of a finite amount of Feynman diagrams. This complies with unitarity order by order 
in the expansion, but fails to fulfill exact unitarity of the scattering amplitude. Thus, some non- 
perturbative resummation should be supplemented to incorporate exact unitarity and, hopefully, 
to accommodate resonances. Regarding this second limitation, several unitarization methods 
have been suggested in the literature and previously used to describe the meson-baryon dynamics: 
Inverse Amplitude Method (IAM) [12], or a somehow modified IAM to account for the large 
baryon masses [13], dispersion relations [14, 15, 16], Lippmann-Schwinger Equation (LSE) and 
Bethe-Salpeter Equation (BSE) [17, 18, 19, 20, 21, 22, 23, 24, 25, 26]. 

In this work, we will study the s— wave meson-baryon scattering for the strangeness S — — 1 and 
isospin 1 = sector in a Bethe-Salpeter coupled channel formalism incorporating CS. This reaction 
provides a good example of the need of unitarization methods; the recent work of Ref. [27] shows 
that HBChPT to one loop fails completely already at threshold. The KN scattering length turns 
out to have a real part about the same size but with opposite sign and half the imaginary expected 
experimentally [28], due to the the nearby subthreshold A(1405)-resonance. It is important to 
realize from the very beginning that, for this particular resonance, the unitarization program 
implies at least considering two channels, namely 7rS and KN. Thus, one deals with a coupled 
channel problem. Unfortunately, there are no one-loop ChPT calculations incorporating this 
coupled channel physics. Much of the discussion below and the rest of the paper reflects this lack 
of knowledge on the general structure of this scattering amplitude within ChPT, particularly on 
the number of undetermined parameters. Of course, not all counterterms are independent, as 
demanded by crossing and chiral symmetries. Actually, the best practical way to impose these 
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constraints on a unitarized partial wave amplitude is by matching to a ChPT amplitude. The 
reason for this is that the left-cut partial wave analytical structure implied by crossing would 
be automatically incorporated within a one-loop ChPT calculation. This point of view has been 
used in several unitarization approaches for the elastic tcN process [15, 12, 13, 14]. There, the 
ChPT amplitude is known up to one-loop [5, 6, 7, 9] which in HBChPT corresponds to third 
and fourth order. This is the lowest order approximation incorporating the perturbative unitarity 
correction which is required to match the perturbative amplitude to a unitarized one. As we have 
already mentioned, the coupled channel ChPT one-loop amplitude for meson-baryon scattering is 
not known, and thus the matching is not possible. 

The first study of the strangeness S — — 1 and isospin 1 = meson-baryon channel incor- 
porating CS and coupled channel unitarization was carried out in Ref. [18, 19], though some 
phenomenological form factors were employed. More recently, this channel has been studied in 
Refs. [21] and [22], where a three-momentum cut-off is used to renormalize the LSE and the off- 
shell behavior is partially taken into account. In principle, the minimal renormalization procedure 
used in Refs. [21, 22] is acceptable but may turn out to be too restrictive in practice. If, instead, 
one takes advantage of the flexibility allowed by the renormalization of an Effective Field Theory 
(EFT), there is a chance of improving the description from threshold up to higher energies. In 
practical terms this means increasing the number of counter terms which one has to add to make 
finite the amplitudes. Along the lines of Ref. [25] , we will adopt this viewpoint below by taking 
fully into account the off-shellness suggested by the BSE using the tree level amplitude as the low- 
est order approximation to the potential. In the approach of Ref. [25] there are three parameters 
for each channel. In the absence of a one-loop ChPT coupled channel amplitude, we content our- 
selves with matching to the tree level one 1 . One could also match the coupled channel unitarized 
amplitude to HBChPT, much below the inelastic thresholds. This procedure was investigated in 
Ref. [25] for the nN system, but did not introduced any powerful constraints. 

Despite of the great success of the model of Ref. [21] in describing the data around the antikaon- 
nucleon threshold, including the features of the Soi — A(1405) resonance, its predictions for higher 
energies [22] do not work so well and clear discrepancies with data appear. Indeed, the limitations 
of the model of Ref. [21] did already appear in the strangeness 5 = and isospin 1=1/2 meson- 
baryon sector, where the model is able to describe data only in a more or less narrow energy 
window around the iV(1535) resonance [23, 24]. However, the model previously developed by two 
of us for the latter channel in Ref. [25] describes data in a wider energy region, ranging from nN 
threshold up to almost a Center of Mass (CM) meson-baryon energy of y/s = 2 GeV, including 
the features of a second resonance (7V(1650)). Motivated by these encouraging results we extend 
in this work the model of Ref. [25] to the strangeness S — — 1 and isospin 1 = meson-baryon 
channel. Like in the S = sector, taking into account the off-shellness of the BSE generates 
a rich structure of unknown constants which allow for a better description of the amplitudes. 
Although the generation of more undetermined constants may appear a less predictive approach 
than putting a cut-off (one single parameter) to regularize the divergent integrals, it reflects the 
real state of the art of our lack of knowledge on the underlying QCD dynamics. The number 
of adjustable Low Energy Constants (LEC's) should not be smaller than those allowed by the 
symmetry; this is the only way both to falsify all possible theories embodying the same symmetry 
principles and to widen the energy interval which is being described. Limiting such a rich structure 
allowed by CS results in a poor description of experimental data. However, a possible redundancy 

1 This point was discussed at length in Ref. [29] in the context of elastic irir scattering; from there one easily 
realizes that the renormalization scheme pursued in Ref. [21] would lead to a scenario where only one of the four 
SU(2)— Gasser-Leutwyler h,l2,h and I4 parameters is independent. 
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of parameters is obviously undesirable but may be detected through statistical considerations (see 
below). The number of LEC's is controlled to any order of the calculation by crossing symmetry. 
In a unitarized approach, the only way to avoid this parameter redundancy is to match the 
unitarized amplitude to the one obtained from a Lagrangian formalism 2 . There is no standard 
one loop ChPT calculation for the meson-baryon reaction with open channels to compare with. 
Some results there exist within the HBChPT scheme, up to order 0(1/M 3 f 2 , 1/Mf*) - being M 
a typical baryon mass - but only involving pions and nucleons [7] . An indirect way to detect such 
a parameter redundancy might be through a fit to experimental data if the errors and correlations 
in some parameters turn out to be very large. We will adopt this point of view in this work, and 
will show that indeed correlations take place effectively reducing the total number of independent 
parameters. 

In this paper four coupled channels have been considered: 7rS, KN, r/A and KE and taken 
into account 577(3)— breaking symmetry effects but neglected the considerably smaller isospin 
violation ones. We also neglect three body channels, mainly for technical reasons. According to 
Ref. [33] the A(1670) resonance decays into the Y,*irN three body final state with a small (.08=F.06) 
branching ratio. 

Preliminary results of the present work can be found in [30, 31, 32]. 

The paper is organized as follows: In Sect, 2, a summary of the theoretical framework used is 
described. Further details can be found in Ref. [25]. Numerical results for amplitudes and cross 
sections and details of the fitting procedure are given in Sect. 3. A special attention is paid both 
to the analytical properties, in the complex plane, of the found t matrix and to the statistical 
correlations between the fitted LEC's, being the Appendices A and B devoted to these issues as 
well. Finally, in Sect. 4 we outline the conclusions of our work. 

2 Theoretical framework 

The coupled channel scattering amplitude for the baryon-meson process in the isospin channel 
/ = 

B(M A , P-k, s A ) + M(m A , k) -> B(M B , P - k', s B ) + M(m B , k') (1) 

with baryon (meson) masses Ma and M B ( itia and m B ) and spin indices (helicity, covariant spin, 
etc..) sa, s b , is given by 

T P [B{k', s B } <- A{k, s A }] = u B {P - k', s B )t P (k, k')u A {P - k, s A ) (2) 

Here, u A {P — k, s A ) and u B (P — k', s B ) are baryon Dirac spinors 3 for the ingoing and outgoing 
baryons respectively, P is the conserved total four momentum and tp(k,k') is a matrix in the 
Dirac and coupled channel spaces. On the mass shell the parity and Lorentz invariant amplitude 
tp can be written as: 

t P (k, k') | on _ shc n = t^s, t)P + t 2 (s, t) (3) 
with s = P 2 = f 2 , t — (k — k') 2 and t\ and t 2 matrices in the coupled channel space. 

2 See discussion in Appendix D of Ref. [25] 
3 We use the normalization uu = 2M. 
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Figure 1: Diagrams summed by the Bethc Salpetcr equation. Kinematics defined in the main text. 
In terms of the matrices t\ and t 2 defined in Eq. (3), the s— wave coupled-channel matrix, 

i 

fo(s) [//], is given by: 
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(5) 



with X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz and Ea,b the baryon CM energies. The phase 

of the matrix Tp is such that the relation between the diagonal elements (A = B) in the coupled 
i 

channel space of / 2 (s) and the inelasticities (77) and phase-shifts (5) is the usual one, 
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(6) 



Further details on normalizations and definitions of the amplitudes can be seen in Section IIB 
of Ref. [25]. 

To compute the amplitude tp we solve the BSE (see Fig. (1)) 



t P (k, k') = v P (k, k') + ij Sjktpb, k')A(q)S(P - q)v P (k, q) 



(2k) 
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where tp(k, k') is the scattering amplitude defined in Eq. (2), vp(k, k') the two particle irreducible 
Green's function (or potential ), and S(P — q) and A(q) the baryon and meson exact propagators 
respectively. The above equation turns out to be a matrix one, both in the coupled channel and 
Dirac spaces. For any choice of the potential vp(k, k'), the resulting scattering amplitude tp(k, k') 
fulfills the coupled channel unitarity condition, discussed in Eq. (21) of Ref. [25]. The BSE requires 
some input potential and baryon and meson propagators to be solved. To compute the lowest 
order of the BSE-based expansion [29] is enough to approximate the iterated potential by the chiral 
expansion lowest order meson-baryon amplitudes in the desired strangeness and isospin channel, 
and the intermediate particle propagators by the free ones (which are diagonal in the coupled 
channel space). From the meson-baryon chiral Lagrangian [2] (see Sect. IIA of Ref. [25]), one gets 
at lowest order for the potential: 



v P (k,k') = t { p\k, k! 
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with D the coupled-channel matrix, 
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(9) 



given in the isospin basis and the same phase conventions as in Ref. [25]. 

While amplitudes follow the chiral symmetry breaking pattern from the effective Lagrangian to 
a good approximation, it is well known that physical mass splittings have an important influence 
when calculating the reaction phase space. Indeed, the correct location of reaction thresholds 
requires taking physical masses for the corresponding reaction channels. We have taken into 
account this effect in our numerical calculation. We also incorporate explicit CS effects to the 
weak meson decay constants and different numerical values for f n , fx and f v can be used. This 
can be easily accomplished through the prescription 



D/f - f-'Df- 



f = Bi&g(f K ,f n ,fr,,f K ) 



(10) 



The solution of the BSE with the kernel specified above can be found in Ref. [25]. It turns out 
that the functions ti and t2, defined in Eq. (3), do not depend on the Mandelstam variable t, and 
thus the dynamics is governed by the matrix function t(s) (see Eq. (4)) 

t(s) = yfit 1 (s)+t 2 (s), (11) 

which is given in Eq. (34) of Ref. [25]. There, the renormalization of the obtained amplitudes 
is studied at length. As a result of the renormalization procedure, and besides the physical and 
weak meson decay constants, a total amount of 12 
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undetermined LEC's appear. We fit these constants to data, as we will see in the next section, 
and from them we define the three following diagonal matrices: 
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which appears in the solution of the BSE. We have denoted the meson-baryon low energy constants 
J (s = (mi + Mj) 2 ), i = K,n,ri,K and j = N, £, A, S of Eq. (A8) of Ref. [25] as J -. 

3 Numerical results 

Throughout the paper we will use the following numerical values for masses and weak decay 
constants of the pseudoscalar mesons (all in MeV), 

m K = nix = 493.68 m n = 139.57 m v = 547.3 

M p = 938.27 M s = 1189.37 M\ — 1115.68 M s = 1318.0 

U = fv = Sk = 1.15 x 93.0 (14) 

where for the weak meson decay constants we take for all channels an averaged value 4 . This 
selection of the coupling constants does not omit the essential features of the meson-baryon system 
since a recent work [34] shows that the observed SU(3) breaking in meson-baryon scatterings 
cannot be explained by the present SU(3) breaking interactions and that the essential physics of 
the resonances seems to lie in the substraction constants. 

3.1 Fitting procedure 

We perform a x 2 ~ fit, with 12 free parameters, to the following set of experimental data and 
conditions: 

• 5 , oi(5 , 2T2j) KN — > KN and KN — > 7rS scattering amplitudes (real and imaginary parts) [35], 
in the CM energy range of 1480 MeV < y/s < 1750 MeV ~ ~ 

In this CM energy region, there are a total number of 56 data points (28 real and 28 
imaginary parts) for each channel. The normalization used in Ref. [35] is different of that 
used here and their amplitudes, T$° 77 , are related to ours by: 

(15) 

where sig(i,j) is +1 for the elastic channel and —1 for the KN — > 7rS one. On the other 
hand, and because in Ref. [35] errors are not provided, we have taken for those amplitudes 
errors given by 

5T?° 77 = ^(O.UT^y + Om 2 (16) 
in the spirit of those used in Ref. [36]. 

• Sqi — ttS mass spectrum [37], 1330 MeV < y/s < 1440 MeV: In this CM energy region, there 
are a total of thirteen 10 MeV bins and the experimental data are given in arbitrary units. 

4 This makes the comparison with Ref. [21, 22] more straightforward. One could also take the physical values 
for these decay constants, but part of the effect can be absorbed into a redefinition of the parameters in Eq. (13). 
See also Eq. (34) in Ref. [25]. 



7f° 77 = sig(*, j)N fo 2 (s) 
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To compare with data, taking into account the experimental acceptance of 10 MeV, we 
compute: 

A(T /•M 7rE (i)+5MeV 
A[M 7rS (i)] J M„ s {i)-5UeV 

where C is an arbitrary global normalization factor 5 and % denotes the bin with central CM 
energy M^(i). Hence, there are only 12 independent data points. Finally, we take the error 
of the number of counts, iVj, of the bin % to be \. f o\\[W i as in Ref. [38]. 

• The K~p -> r]A total cross section of Ref. [39], 1662 MeV < y/s < 1684 MeV : We use the Crys- 
tal Ball Collaboration precise new total cross-section measurements (a total of 17 data points 
compiled in Table I of Ref. [39]) for the near-threshold reaction K~p — > r)A, which is domi- 
nated by the A(1670) resonance. We assume, as in Ref. [39], that the p— and higher wave 
contributions do not contribute to the total cross-section. 



fHs = x 2 ) 



2^2 



\k 2 (s = x 2 )\ x 2 dx 



(17) 



Finally, we define the x 2 5 which is minimized, as 

1 N -, n a / (a)tfc (a)\ 2 

^ a =i n a j=1 \ a) J 

where N = 4 stands for the four sets of data used and discussed above 6 and x^ th denotes 
our model result for the data point Xj*\ Finally n a takes the values 56, 56, 12 and 17, and 
N t ot = E«=i n a is the total number of data points. With such a definition, acceptable best fits 
should provide values of xV^tot around one. 

Though we have considered four coupled channels, three-body channels, for instance the -kttH 
one, are not explicitly considered, as it has been also assumed previously in Refs. [22] and [39]. 



3.2 Results of the best x 2 fit 

The model presented up to now has initially twelve free parameters (Eq. (12)), which have to be 
determined from data. This is a cumbersome task because there are many mathematical minima 
which are not physically admissible. For instance, in some cases one finds fits to data with spurious 
poles in the first Riemann Sheet which strongly influence the scattering region and hence violate 
causality. Any fit embodying these singularities is physically inadmissible and should be rejected. 
This is an important issue, which should always be considered in any analysis. We will further 
elaborate on this point in the Appendix B. 

Furthermore, even if one is reasonably convinced that a physically acceptable minimum has 
been found, there are strong correlations between the fitted parameters, which have to be carefully 
evaluated and, if possible, understood. Our best results come from a minimum for which the pairs 
( Jj, A B .) with % = KN, 7rS, r/A, KE and B, L = N, S, A, S are totally correlated (correlation factors 
bigger than 0.99) which leads to an almost singular correlation matrix reflecting the fact that 
there exists, in very good approximation, a linear relation between the Jj and A B . parameters 7 . 

5 We fix it by setting the area of our theoretical spectrum, J2i a[m ' to tne tota l number of experimental 
counts J2i N i- 

6 From the first item above and to define the x 2 , we consider two separated sets: K N — > KN and KN — > 7rS. 
7 In Ref. [25] the static limit (infinitely heavy baryons) is discussed, and it is shown (Eq. (D9)) that there exists 
a linear relation between Jj and A^, if, as is the case here, A TOi , i = K, ir, r/, K is small when compared to . 
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Thus, we have fixed Ab x to some specific values in the neighborhood of the minimum, given in 
Appendix A, and have studied the correlation matrix for the remaining eight parameters. Yet, 
we find a strong correlation (0.99) between Jke and Ak and we proceed as above, i.e., we fix Ak 
and evaluate the correlation matrix and variances for the remaining seven parameters. Thus, at 
the end of the day we have only seven independent best fit parameters. The best fit parameters, 
their variances and the correlation matrix are compiled in Appendix A. 

In Figs. (2-4) we compare the results of our best fit with the experimental data. The overall 
description is remarkably good, both at low energies and the higher end of the considered energy 
region. Besides, as we will see, the description of the A(1405) and A(1670) features is also quite 
good. Thus, our scheme leads to a much better description of the data than the approach of 
Ref. [22], as it was also the case in the strangeness S = sector ([25] versus [23]). 




1320 1340 1360 1380 1400 1420 1440 1460 1660 1665 1670 1675 1680 1685 1690 

E CM [MeV] E CM [MeV] 



Figure 2: Solid lines: Results of our fit. Experimental data for 7rS — > 7rS and K p — > rjA are from 
Refs. [37] and [39], respectively 




Figure 3: The real (left panel) and imaginary (right panel) parts of the s— wave T— matrix, with normalization 
specified in Eq. (16), for elastic KN — ► KN process in the I — isospin channel as functions of the CM energy. 
The solid line is the result of our fit, and the experimental data are taken from the analysis of Ref. [35] with the 
errors stated in the main text. 
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Figure 4: Same as in Fig. 3 for the inelastic channel KN — > 7rS. 



For the elastic KN — > KN scattering length we get 



a KN 



f~ 2 (s = (m K + M, 



KN^KN 



= (-1.20 ± 0.09 + i 1.29 ± 0.09) fm 



(19) 



where the error is statistical and it has been obtained from the covariance matrix given in the 
Appendix A, taking into account the existing statistical correlations, through a Monte-Carlo 
simulation. This value should be compared both to the experimental one (—1.71 + i 0.68) fm of 
Ref. [28] and to the LSE approach of Ref. [21] (—2.24 + i 1.94) fm. Unfortunately, the previous 
works do not provide error estimates, so one cannot decide on the compatibility of results. 



3.3 Second Riemann Sheet: poles and resonances. 

In this section we are interested in describing masses and widths of the Sqi — resonances in the 
S — — 1 channel. Since causality imposes the absence of poles in the t(s) matrix in the physical 
Sheet [40], one should search for complex poles in unphysical ones. Among all of them, those 
closest to the physical Sheet and hence to the scattering line are the most relevant ones. We 
define the second Riemann Sheet in the relevant fourth quadrant as that which is obtained by 
continuity across each of the four unitarity cuts (see a detailed discussion in a similar context 
in Ref. [25]). Physical resonances appear in the second Riemann Sheet of all matrix elements of 
t(s), defined in Eq. (11), in the coupled channel space, differing only on the value of the residue 
at the pole. The residue determines the coupling of the resonances to the given channel. In 
Fig. 5 we show the absolute value of the i]A — > rjA element of the t matrix. We choose this 
channel because all found poles have a sizeable coupling to it. Both the fourth quadrant of the 
second Riemann Sheet and the first quadrant of the first (physical) Riemann Sheet are shown. 
The physical scattering takes place in the scattering line in the plot (upper lip of unitarity cut of 
the first Riemann Sheet). We find three poles in the second Riemann Sheet which positions are 
(s = M 2 R -iM R T R ): 
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Figure 5: Modulus of the rjA — > ?7A element of the scattering amplitude t(s) [fm], defined in Eq. (11), analytically 
extended to the first and fourth quadrants of the s— complex plane. The solid line is the scattering line, s = 
x + i + , from the first threshold, (m^ + Ms) 2 , on. The first (second) Riemann Sheet is depicted in the first 

(fourth) quadrant of the s complex plane. Three poles appear in the second Riemann Sheet, which are connected 
with the A(1405) and A(1670) resonances, see discussion in the main text. Besides unphysical poles show up in 
the physical Sheet out of the real axis, but they do not influence the scattering line as can be seen in the plot. 
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First Pole : M R = 1368 ±12 T R = 250 ± 23 (20) 

Second Pole : M R = 1443 ±3 T R = 50 ± 7 (21) 

Third Pole : M R = 1677.5 ± 0.8 T R = 29.2 ± 1.4 (22) 

where all units are given in MeV and errors have been transported from those in the best fit 
parameters (Eq. (45)), taking into account the existing statistical correlations through a Monte- 
Carlo simulation. 

These poles are related to the two S i resonances A(1405) and A(1670) which appear up to 
this range of energy in the PDG (Ref. [41]). The third pole above can be clearly identified to the 
A(1670) which is located at 

A(1670) : M R = 1670 ±10 T R = 35±$o Ref - [ 41 ] 

M R = 1673 ±2 = 23 ± 6 Ref. [33] (23) 

where again units are in MeV. The agreement of our predictions and the experimental data is 
satisfactory and better than the previous theoretical LSE approach of Ref. [22]. Let us look at 
the A(1405) resonance, which nature is under much discussion [42, 43]. Following the PDG it is 
placed at (in MeV) 

A(1405) : Mr = 1406.5 ± 4.0 T R = 50 ± 2 Ref. [41] (24) 

Our amplitudes have two poles in the region of 1400 MeV, Eqs. (20) and (21). The features 
of the second one are in agreement with the previous results of Refs. [21, 22] and though the 
width compares well with the experiment, the mass is shifted to higher values. Besides, we should 
note that the pole quoted in Eq. (20) is very broad and can not be identified with any of the 
experimentally established resonances. This pole is also present in the LSE model of Refs. [21, 22], 
as it was pointed out in Ref. [44], though the mass position there is similar (M R = 1390 MeV), the 
width is about a factor two narrower (r^ = 132 MeV) than ours. Our understanding is that this 
broad resonance does not influence strongly the scattering line. However, the 7rS mass spectrum 
peaks around 1405 MeV in the experimental data and also in our approach as can be seen in 
Fig. 2. This is a clear indication of a sizeable non resonant contribution on top of our 1443 MeV 
pole. 

On the other hand, there are unphysical poles in the physical (first) Riemann Sheet. These 
unphysical poles appear because we have truncated the iterated potential to solve the BSE. The 
two of them closer to the scattering line are located at (s = M 2 + iMT) with M pa 1166, T pa 
±200 MeV and M pa 1616, r pa 631 MeV. The tails of both poles can be seen in Fig. 5 and they 
do not influence the scattering line. In Appendix B, we will show the results from a fit which, at 
first sight, are even in a better agreement with the experimental data (Subsect. 3.1) than those 
presented up to now. However this apparent improvement is achieved because the unphysical poles 
get closer to the real s axis and they affect, in a substantial manner, the scattering amplitudes. 
Hence, we discard this minimum, and we would like to note that it is important to observe the 
positions and influence of the unphysical poles when deciding the goodness of a phenomenological 
description of data. 
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Finally, we have also analyzed the nature of resonances on the light of the well known Breit- 
Wigner parameterization for coupled channels and real s (See e.g. Ref. [45] and references therein), 



X r n e m+Sj) M /pBWpBW 

= ^ - 1 + !_ r (25) 



2i Pi L J v ^"p- [s-M 2 R + yMrTr 



for which the background is assumed to be diagonal in coupled channel space and the relative 
phase of the resonance to the background and the sum partial decay widths, X^rf w = T R , 
are chosen in such a way that t BW (s) exactly fulfills unitarity on the real axis. Here, p« is a 
kinematic factor defined by the second line of Eq. (31) below. The branching ratio is then defined 
as Bf w = rf w /T R . Subtracting the resonance contribution, Eq. (26), to the total amplitude we 
have found that for our A(1670) the background is not a diagonal matrix, since for our ^-matrix 
we get 2 Ei<i \Uj ~ ^ BW | 2 ~ Ei \Ui ~ tT W \ 2 f o r s -> M%, with t^ BW the second term in Eq. (25). 
In addition, the BW parameterization suggests a relation between the residue at the pole and 
the imaginary part of the pole. This relation is only true in the sharp resonance approximation, 
pBW << p. w j^] 1 p. th e CM momentum of the decaying state. We have also checked that for our 
problem this is not the case. Actually, with such a definition we find that J2i Tf W ~ 0.81^ for the 
A(1670). This is a simple consequence of the incorrect assumption made by Eq. (25). 



3.4 Branching ratios and couplings of the resonances to different final 
states. 

Before going further we would like to make some critical remarks regarding both the comparison 
between "theory" and "experiment". Our BSE solution has a very specific energy dependence 
which, as we have seen in Subsect. 3.2, is able to fit numerically experimental data or rather a 
partial wave analysis with a given energy dependence. Obviously, both functional forms are not 
identical, and it is also fair to say that both incorporate their own biases. There is no reason 
to expect that they are numerically alike also in the complex plane 8 . Under these conditions, 
some parameters, like branching ratios, have a different meaning, since the extrapolation of the 
resonant contribution to the real s-axis is ambiguous. Actually, the ambiguity is enhanced as the 
resonance becomes wider and as a consequence the definition of a branching ratio becomes model 
dependent. We explain below our definition of branching ratios and how they are extracted from 
our amplitude. 

Let us consider s R = M R — i M R Y R a pole in the second Riemann Sheet of the coupled channel 
scattering matrix t(s). Then, around the pole, it can be approximated by 

[t(s)] » 2M R ^^, (26) 
J s- s R 

where g^ is the residue matrix. Since t is a complex symmetric matrix (due to time reversal 
invariance), g is also complex symmetric and its rank is one to ensure that s = s R is a pole of 
order one of the det[t(s)]. In this way a non-degenerate resonant state is being described 9 . Under 

8 A good example of this fact is provided by our best fit results of Subsect. 3.2 and the physically inadmissible 
results of Appendix B; they look very much the same on the scattering line although the analytic structure is 
rather different. 

9 This can be seen as follows. Using a matrix notation, the BSE reads t(s) = V + VGo(s)t(s) and it is solved by 
t(s) = V(l — Gq(s)V)^ 1 , with the obvious identifications for V and Go. A pole, at s = sr, in det[t(s)] is produced 
by a zero of det[l — Go(s)V]. This last condition ensures that the homogeneous (quasi)bound state Bethe-Salpeter 
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these conditions, turns out to be factorizable 10 , 

9ij = 9%9j- (28) 

The above matrix, g^, has only one non-zero eigenvalue, g\ + g\ + g\ + g\, with g { the associated 
eigenvector. The vector ^ determines the coupling of the resonance to the different final states, 
which are well and unambiguously defined even if the corresponding channels are closed in the 
decay of the resonance. In Table 1 we give the complex vectors gi for the three resonances described 
in Subsect. 3.3. Unfortunately, the PDG does not provide this kind of information and, instead, 
branching ratios are given. To extract meaningful branching ratios from our calculation, we have 
to extrapolate the resonant contribution of the scattering amplitude to the s real axis, which is the 
only experimentally accessible. In addition, the picture of a resonance as a quantum mechanical 
decaying state, requires a probabilistic description. Thus, we isolate the resonant contribution to 
the S-matrix 11 for s = M R 

S—\s = M 2 R ) = -2i2M R ^-§^L- 

M R - S R 



Pl (s) = e( s -4)^(v^T^ + ^) 



Pl (s = M 2 R ) (31) 



with s\ h the threshold of the baryon-meson channel i. This definition embodies a sensible kinematic 
suppression compatible with Cutkosky's rules and the s-wave nature of the resonance. Defining 




(32) 



we find 



S~\s = Ml) = -2i M R V R J ibj (33) 

M R - S R 

and taking into account that the matrix bij = bfij has rank 1 and that bi (or any vector proportional 
to it) is the only eigenvector of 5' resonant with a non zero eigenvalue, the resonant state at s = M R 

equation 

(-Go(s)- 1 + V)V = (27) 

has a non trivial solution for s = s^. Indeed, all solutions, of the above equation are linear combination of 
the null eigenvectors of the (1 — Go(sr)V) matrix and describe the dynamics of the existing states (resonances) at 
s = sr. For a non-degenerate resonance (which we have checked it is indeed the case), the zero eigenspace should 
have dimension one. Whence for the case of four coupled channels, the rank of the matrix (1 — Go(sr)V) is three 
and therefore det[l — Go{s)V] should have a single zero at s — sr. 

10 The symmetric complex matrix g can be diagonalized by a complex orthogonal transformation, U, as g = 
U T d U, being only one element of the diagonal matrix d different of zero. If we take this element to be the dn, 
we have = {U T )ndi 1 Ui = dnUuUi-j . 

11 The S— matrix is related to the t matrix, in our convention, by: 



S l3 (s) = 5ij - 2i v^) tij (s) yj Pi (s) (29) 
and probability conservation (S^S = SS* 1 = 1) holds since t fulfills coupled channel unitarity: 

tt.(s) - = 2i Pk (s) t kj (s). (30) 
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9kn 


9t 


-s 


9 V A 




9k~ 


Resonance 


\g\ 


\g\ 





\g\ 





\g\ 





M R = 1368 MeV 


3.9(1) -0.59(5) 


3.65(8) 


-0.73(3) 


1.7(2) 


3.0(2) 


0.29(7) 


1.14(13) 


M R = 1443 MeV 


3.3(2) 0.72(7) 


2.14(16) 


1.10(8) 


2.2(1) 


-2.66(3) 


0.23(1) 


-0.008(57) 


Mr = 1677.5 MeV 


0.39(2) -1.29(4) 


0.20(1) 


0.77(5) 


1.22(3) 


2.69(2) 


1.64(1) 


-0.13(1) 



Table 1: Dimensionless complex couplings, gi = \gi\ e 1 ^, defined in Eqs. (26-28), for all channels i = 
KN, 7rE, r/ A, K3, and for the three poles (resonances) quoted in Eqs. (20-22). The phases are in radian 
units. Errors are purely statistical and affect the last significant digit. 



will be given by 

\R)<xJ2bi\i) (34) 

i 

where \i) stands for the meson-baryon states used to build the coupled channel space. Finally, the 
branching ratio B^ will be given by the probability of finding \R) in the state \i) 

B ' = (35) 

which by definition fulfills Y,i Bi = 1. The partial decay width may then be denned as Tj = B^Yr, 
and obviously J2i Tj = F R . For the A(1670) we obtain the following branching ratios with the 
above prescription: 

B RN = 0.24 ±0.01, B^ = 0.08 ± 0.01, B vA = 0.68 ±0.01 (36) 

The last two values are not in agreement with the values quoted in the PDG [41] ( B^ N = 
0.25 ±0.05, B^ = 0.40 ±0.15, B vA = 0.17 ± 0.07 ) and in Ref. [33] ( B kN = 0.37 ± 0.07, S wE = 
0.39 ± 0.08, B vA = 0.16 ± 0.06, Brians) = 0.08 ± 0.06 ). 

To finish this subsection, we would like to point out that in the present context the concept of 
branching ratio is subtle and it might be ambiguous, both from the theoretical and experimental 
sides. From the experimental point of view the difficulty arises from the impossibility of preparing 
a pure short-lived resonant state strongly coupled to a continuum and therefore the impossibility 
of disentangle events coming from the formation of the resonance from those produced through 
non-resonant processes. From the theoretical point of view the ambiguity comes when denning 
^resonant ^ s _ j^2^ j n (31). For instance at the pole s = s R one could have 

A(g) 9ij Pjjs) 
l (3i(s R ) s-s R l3j(s R y 

instead of the expression assumed in Eq. (26), being f3i(s) an arbitrary, analytical around sr, 
complex function. In these circumstances the matrix S*j Sona,nt (s = M R ) would be different from 

that given in Eq. (31) by a factor ^gg and one would get a new vector b; L which in terms 

of the vector bi, defined in Eq. (32), reads: 

leading to, in principle, different branching ratios. The trouble comes from the extrapolation from 
sr to the real axis which is not unique. 



« 2Mr^^-^-^^, (37) 
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Figure 6: J = meson-baryon s— wave cross sections for different channels. 

The usual assumption is that the fii(s) functions are smooth and they do not change much 
from s = sr to s = M|, and more important, that the change does not depend significantly on 
the channel i. However, this is not always true, for example, if one were dealing with a p— wave 
resonance the function f3i(s) would at least include the CM meson baryon momentum. If there is 
a channel which gets open close and above to s = Mj,, then the CM momentum would lead to a 
suppression of the branching ratio to this channel. 

3.5 Predictions for other processes 

In Fig. 6 we show some of our predictions for s— wave 1 = cross sections for some elastic and 
inelastic channels. For most of them there are no data. The effect of the A (1405) resonance is 
clearly visible in the S7r — > S7r and S7r — * KN and KN — > KN cross sections. On the other hand, 
the elastic 77A cross section takes a very large value at threshold, which corresponds to a typical 
low energy resonance behavior triggered by the A(1670) resonance. This is in contrast to any 
expectation based in the Born approximation, since the corresponding potential in this channel 
vanishes (Eq. (9)). Our estimates for 7rS and r/A scattering lengths (defined for the elastic channels 
similarly as in Eq. (19)) are 

a^ s = 1.10 ± 0.06 fm 

a vh = (0.50 ± 0.05) + i (0.27 ±0.01) fm (39) 

respectively. 
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3.6 Heavy Baryon Expansion at Threshold 

As we have already mentioned, the only calculation within HBChPT in the S — — 1 sector is that 
of Ref. [27], where the s-wave scattering lengths in both isospin channels 1 = and 1 = 1 are 
computed. It is found that HBChPT to one loop fails completely in the 1 = channel due to the 
strong influence of the subthreshold A(1450) resonance. We think that it is of interest to analyze 
this problem within the context of our unitarization approach. 

As was discussed previously by two of us [25] the condition for the coupled channel amplitude 
to have a well defined static limit, M — > oo, is that the combination 



CrhM — - 



1 

m 



^4m + ^ A ™- A m) 



(40) 



goes to some definite finite value. The parameters J^m an d A m and A^ are listed in Eq. (13) 
and their best fit result is presented in Eq. (45) of Appendix A. Using Eq. (45) one can estimate 
the combinations in units of the relevant meson mass 

A w -A s 



M N J RN + 
M A J„ A + 
M S J KS + 



4M S 
A R -A N 



AM, 



N 



4M A 
A K -A S 
4M= 



+0.24™* 
-0.09™^ 

+0.08771,, 

-0.64mA- 



(41) 



showing that they are not unnaturally large. Following the discussion of Ref. [25], the heavy baryon 
expansion is done in the standard way, and taking the leading heavy baryon approximation of the 
parameters, 



J. 



mM 
Am 

Am 



= f 



rhM 



A° 



1 + 0| F 



(42) 



We get for the KN s-wave scattering length in the 1 = channel the following expressions 
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(43) 



16 



This expression can be mapped into the HBChPT result of Ref. [27], since the trascendental 
function dependence is exactly the same. This should be so because they build the perturbative 
unitarity correction in HBChPT. Thus, one could identify a linear combination of the leading 
order approximation of our constants, with another linear combination of HBChPT constants. 
On the other hand, if we assume the values of the best fit parameters, Eq. (45), for and 
and A°~ we obtain the following numerical estimate 

tt£ KN r?A 

a RN = L02- 053 + 028 + 1.00 + i 0.41 -L25 + 317 + . . . = 3.69+ i 0.41 f m + . . . (44) 

l// 2 V/ 2 M 1// 2 M2 1//4 

where the contributions are separated according to the order in the chiral expansion and also to 
the corresponding intermediate state. As we see, large cancellations at higher orders must take 
place to obtain, after summing the whole series, our result in Eq. (19). Note also that the real 
part is about twice and with opposite sign as compared to the experimental result [28]. A similar 
situation occurs in HBChPT [27]; the real part of the scattering amplitude has opposite sign 
although similar magnitude to the experimental number. Likewise, large cancellations have also 
been noted, indicating a bad convergence rate in HBChPT. The fact that our calculation, Eq. (44), 
gives a larger magnitude for Re a^ N than in HBChPT reflects, in addition, a bad convergence in 
the expansion (43). This situation has also been described in the coupled channel case S = 
sector, [25] and seems a common feature of unitarization methods [12, 13, 14]. 

4 Conclusions 

In this paper we have extended the Bethe-Salpeter formalism developed in Ref. [25] to study 
s— wave and 1 = meson-baryon scattering up to 1.75 GeV in the strangeness S — — 1 sector. We 
work on a four dimensional two body channel space and the kernel of the BSE takes into account 
CS constraints as deduced from the corresponding effective Lagrangian. The obtained t matrix 
manifestly complies with coupled channel unitarity and the undetermined low-energy constants 
of the model have been fitted to data. The available direct experimental information is limited 
to the 7rS — > 7rS mass spectrum, and the K~~p — > Arj total cross section, for which errors are 
provided, and to KN — > KN and KN — > 7rS scattering amplitudes of a partial wave analysis, for 
which errors are guessed. Taking this into account, the agreement with experiment is satisfactory. 
Besides, some predictions for other cross sections, not yet measured, have been also given. A 
careful and detailed statistical study has been carried out, showing that only seven parameters 
(LEC's) out of the starting twelve are really independent. Thus, though our model has some more 
free parameters than those required in Ref. [22], the description of data achieved in our approach 
is superior to that of Ref. [22] . A similar situation was already found also in the strangeness 5 = 
sector [25]. According to previous experience, the reduction of parameters is partly due to the 
constraint of a well defined heavy baryon limit [25]. Likewise, crossing symmetry is expected to 
shed more light on the number of independent LEC's. As we have argued, matching to HBChPT 
calculations in the S — — 1 and 1 = sector would be the ideal way to map our LEC's into those 
stemming from an effective chiral Lagrangian, but it is already known that the chiral expansion 
fails [27] to reproduce the KN scattering length, due to the influence of the nearby subthreshold 
A(1405) resonance. Possibly this could be overcome by properly accounting for the singularity 
structure as suggested in [15]. We believe these points deserve a deeper investigation. 

We have undertaken a careful discussion on the analytical structure of the scattering matrix 
amplitude in the complex s— plane, which becomes mandatory in order to extract the features 
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of the Sqi resonances. We have searched for poles in the second Riemann Sheet and compared 
masses, widths and branching ratios to data. The agreement is also quite satisfactory. In the 
resonance region our unitary amplitude cannot be analyzed as a Breit-Wigner resonance due to 
a sizeably non-diagonal background in coupled channel space. This, in particular, prevents from 
a simple interpretation of branching ratios. Although residues at the resonance poles are well 
and unambiguously defined, the definition of branching ratios requires special considerations and 
provisos, due to an ambiguous extrapolation of the resonance contribution of the S-matrix from 
the pole to the scattering line. We have also illustrated that looking for a good description of 
experimentally accessible data is not sufficient and that, in some cases, it can be achieved at the 
expense of generating non-physically acceptable poles in the first Riemann Sheet, which influence 
on the scattering region is non-negligible. Thus, any fit to data should be supplemented by 
this additional requirement of not producing spurious singularities numerically relevant for the 
description of scattering processes. 
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A Best fit results 

The best fit (x 2 /N tot = 0.93) parameters are 







Jrn 


= -0.0186 


±0.0010 
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±0.00061 
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±0.00021 








= -0.11936 


±0.00018 
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+ M S ) 2 
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±0.00036 


A A 


= A A / (m v 


+ Ma) 2 


= -0.00262 


±0.00011 



(45) 



with fixed parameters 

A s /(m K + M E ) 2 = -0.0035 

A R /{m R + M N ) 2 = -0.034 

A n /(m n ±M S ) 2 = 0.060 (46) 

A v /(m v + M A ) 2 = 0.049 

A^/(m^±M H ) 2 = -0.26 

as explained in the main text. We assume that the parameters of Eq. (45) are Gaussian correlated, 
this is justified because they come from a \ 2 ~ fit- To make any further statistical analysis of 
quantities derived from the parameters above, the corresponding covariance (v ) and correlation 
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(c) matrices are needed. These matrices are defined as usual 



Vij = 



Cij — 



1 dx 2 " 
2db k dbi , 



J ij 



being any of the seven parameters J's and A's of Eq. (45). The errors, Sbi 
are obtained from the diagonal elements of the covariance matrix ( 5bi = 
estimate for the correlation matrix reads: 



(47) 

quoted in Eq. (45) 
y/vii). Finally our 
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Despite the correlations in the above matrix are at most in modulus of about 0.9, the matrix has 
got an eigenvalue quite close to zero (0.0025), which is a clear indication that one of the parameters 
might still be redundant. 



B Non-physically acceptable fits to data 

In Figs. (7-9) we present the results of a fit to the data, which we will show is not physically 
acceptable. The overall description of data is remarkably good. This fit gives xV-^tot = 0.69 to 
compare with the value of 0.93 of the best fit presented in the main text. However, when 
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Figure 7: Same as in Fig. (2) for "the non physically acceptable fit" described in the Appendix B. 



looking at the t(s) matrix in the s— complex plane (Fig. (10)), one realizes that there exists a 
proliferation of poles, some of them unphysical and others with no experimental counterparts. In 
the first Riemann Sheet we find at least two poles. The first one is located at s = M 2 — iMT with 
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Figure 8: Same as in Fig. (3) for "the non physically acceptable fit" described in the Appendix B. 
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Figure 9: Same as in Fig. 8 for the inelastic channel KN — > 7rS. 



M = 1606 MeV and T = 153 MeV. This pole is close to the s real axis and produces visible effects 
on the scattering line not only for the i^A-^A entry shown in the figure, but also for all % — > j 
channels. Indeed, as we showed in Subsect. 3.3, our preferred fit has also a similar unphysical pole 
but significantly farther (r = 631 MeV) from the real axis, and therefore with a tiny influence on 
the physical scattering. Since causality imposes the absence of poles in the physical Sheet [40], 
the existence of such a pole affecting the scattering line invalidates the description ( Figs. (7-9)) 
presented in this Appendix, despite of its quality. Besides in the first Sheet, there exists a pole 
on the real axis (-^/s = 1307 MeV) and below the first threshold, which would correspond to a 
bound state, stable under strong interactions. Such a state has the same quantum numbers as the 
A(1115) and it should show up in all reactions where the latter one is produced. 

On the other hand, in the second Sheet there exist now four poles. Three of them are similar 
to those presented in Subsect. 3.3, though the one placed around 1370 MeV (Eq. (20)) is now 
almost a factor two narrower (it is located at Mr = 1392 MeV and Tp> = 120 MeV). In addition 
there exists a new resonance Mr = 1343 MeV and Tr = 0.18 MeV which is responsible of the high 
peak at the beginning of the scattering line in Fig. 10 and of the existing bump between 1330 and 
1360 MeV in the 7rS mass spectrum of Fig. 7. As far as we know, there are no other independent 
indications of the existence of this extremely narrow resonance. 

We have presented the results of this " non physically acceptable fit" to stress that, in order 
to be sure of having a good approach to the t matrix of a given physical system, one should not 
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Figure 10: Same as in Fig. (5) for a the non physically acceptable fit to data described in the AppcndixB. 

only look at the t matrix and related observables (cross sections, ...) at the physical scattering 
line (real s values), but one should also study its behavior on the s— complex plane, both in the 
second Riemann Sheet to find the resonances and in the first Riemann Sheet to be sure of avoiding 
pathological behaviors as the one illustrated in Fig. 10. 
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